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Abstract. The prompt emission from gamma-ray bursts (GRBs) still requires a physical explana- 
tion. Studies of time-resolved GRB spectra, observed in the keV-MeV range, show that a hybrid 
model consisting of two components, a photospheric and a non-thermal component, in many cases 
fits bright, single-pulsed bursts as well as, and in some instances even better than, the Band func- 
tion. With an energy coverage from 8 keV up to 300 GeV, GLAST will give us an unprecedented 
opportunity to further investigate the nature of the prompt emission. In particular, it will give us the 
possibility to determine whether a photospheric component is the determining feature of the spec- 
trum or not. Here we present a short study of the ability of GLAST to detect such a photospheric 
component in the sub-MeV range for typical bursts, using simulation tools developed within the 
GLAST science collaboration. 
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1. INTRODUCTION 

The Band function [see 1], a softly broken power-law, describes the keV-MeV spectra 
for a broad range of gamma-ray bursts (GRB) exceptionally well. It is an empirical func- 
tion without any claim to describe the physical processes behind the continuum spectra. 
However, it has properties that in many bursts could be described as the result of non- 
thermal radiation processes, such as synchrotron or inverse Compton scattering. Some 
spectra, even though well-described by the Band function, can not be the result of plain 
optically-thin synchrotron (OTS) processes since the photon index of the sub-peak spec- 
tral slope have a value a larger than —2/3, and hence lies beyond the synchrotron "line 
of death", LOD [see 2, 3]. Synchrotron processes may however still be a viable option 
through realization of synchrotron self-absorption (SSA) in a higher energy domain, 
see e.g. Lloyd and Petrosian [4]. The physical conditions required to produce such hard 
spectral slopes in the 100 keV domain are generally considered to be unreasonable but 
too little is known regarding the generation of magnetic fields and the relevant condi- 
tions in order to completely rule out this possibility. Inverse Compton scattering with 
a self-absorbed synchrotron seed spectrum residing in the optical domain {synchrotron 
self-Compton, SSC) is a possible mechanism for spectra with a larger or equal to 0, see 
Panaitescu and Meszaros [5]. Still another explanation for the hard spectral slopes, sug- 
gested by Medvedev [6], is a combination of jitter radiation and OTS radiation. Baring 
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and Braby [7] argued, however, that since both the latter emission processes need an al- 
most purely non-thermal electron distribution to be able to fit the observed spectra, these 
emission models are difficult to reconcile with the assumed diffusive shock acceleration 
models, which often give rise to a strong contribution of a thermal population. 

Models predicting a photospheric emission component in the prompt spectra of GRBs 
has been suggested by several authors [see 8, 9, 10, 11, 12, 13, 14, 15]. Meszaros 
and Rees [8] argue that the contribution of thermal radiation - originating from the 
expanding fireball when it becomes optically thin to Thomson scattering - could explain 
the hard spectral slopes in the prompt GRB spectra. Ryde shows in [12] that some bright 
CGRO BATSE bursts have sub-peak spectral slopes that fit blackbody radiation very 
well, in agreement with results by e.g. Ghirlanda et al. [16], and also demonstrates in 
[13] that time-resolved BATSE spectra can, as an alternative to the phenomenological 
Band function, be fitted with a hybrid model consisting of a single power-law and a 
blackbody function. The hybrid model gives, for bright single pulsed bursts, just as good 
fit and in many cases even a better fit than the Band function. Hence, by combining a 
photospheric component with a power-law function, that can be interpreted as being a 
result of further energy dissipation in the optically thin part of the outflow, a wide range 
of GRBs observed by BATSE can be explained. 

In this paper, we investigate the GLAST detectability of bursts containing a strong 
photospheric component in the sub-MeV range. A sample of 847 time-resolved spectra 
from 57 bright BATSE bursts were analyzed where both the Band function and the 
hybrid model was imposed on the data. The fits and analysis were performed with 
XSPEC [17]. Three representative bursts were selected for further detailed analysis 
and used as basis for simulations in the expanded energy range of the two instruments 
onboard GLAST: the GLAST Burst Monitor (GBM), covering an energy range from 
8 keV to 30 MeV (hence encompassing energies both below and above the BATSE 
window), and the Large Area Telescope (LAT), including energies from 20 MeV up to 
300 GeV. The three selected bursts, GRB911016, GRB941026 and GRB960530, were 
all single-pulsed, had comparable ;^^-values and "goodness-of-fit" (P-values) for both 
models, and showed a low-energy power-law index (Band function) that lied beyond 
LOD, i.e. with a larger than —2/3. The temporal features associated with each of the 
three bursts were extracted and the resulting data extended from the BATSE energy range 
(20 keV - 2 MeV) into the GLAST domain (8 keV - 300 GeV), hence assuming a strong 
high-energy emission. The high-energy part is generally considered to be characterized 
by the Band function spectral index j8 that for the majority of BATSE bursts has a value 
less than —2. Of the few observed bursts that has high-energy spectra available, some do 
however show a strong high-energy emission [18]. In section 2 we first describe how the 
BATSE data was parameterized and the model that was used for the GLAST simulations. 
Section 3 explains how the simulated GLAST data was produced and section 4 presents 
the results. We conclude this paper with a discussion in section 5. 
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2. MODELING THE THERMAL GRBS 



2.1. The hybrid model 

The hybrid model consists of two components, a photospheric and a non-thermal 
radiation spectrum. These were in our analysis and in the subsequent simulations repre- 
sented by a blackbody spectrum (Planck function) and a single power-law, respectively. 
This model is in XSPEC described by Nhybrid{E) = Ni,b{E) + Npo{E) where: 

E^dE 

N,,(E) = A„-jjj5^ (1) 

N^E) = A,„ (^) ' (2) 

Ai,i, and A po are normalization constants for respective component, kT the temperature in 
keV for the blackbody component and p the photon index of the power-law component. 
This simple two-component model was used when analyzing the BATSE data in the 
20 keV-1 MeV energy range. The left panel in figure 1 shows a hybrid model fit to 
a time-resolved spectrum of GRB911016 covering the time-interval 0.704-1.280 s after 
the trigger. In our simulations for GLAST we extended the energy range into the GLAST 




FIGURE 1. Time-resolved vFy spectra of GRB9 11016 assumed to have a strong high-energy com- 
ponent in the GeV-domain. See also figure 6. Left panel: XSPEC Hybrid model fit in BATSE energy 
band: 0.704 - 1.280 s. Right panel: SBM hybrid model simulation of the same burst extrapolated in the 
GLAST energy band: 1.0- 1.1 s. The solid line describes the hybrid model, the dashed line the blackbody 
component and the dot-dashed hne the broken power-law. The two vertical lines in the plot describe the 
energy range covered by BATSE. 

domain, and therefore used an extended time-dependent hybrid model that included 
a blackbody function that evolves over time and a broken power-law - instead of a 
single power-law component - with a high-energy cut-off. The blackbody component 
that evolves over time is described by: 

E^dE 

mE.t)=A{tU-^j^^^^^ (3) 
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TABLE 1. Power-law photon indices extracted with BATSE and used as input to SBM 



GRB911016 GRB941026 GRB960530 

Time range [s] Index Time range [s] Index Time range [s] Index 



0.000-3.520 1.30±0.05 0.030-6.272 1.57±0.09 0.029- 18.624 1.42±0.07 
3.520-7.105 1.69±0.08 6.272- 14.784 2.12±0.04 



We motivated the broken power-law based on the following two facts: a) OTS radiation 
may according to theory show a break in the energy domain covered by GLAST and b) 
power-law index values in the BATSE data from the hybrid model fits can be interpreted 
as being distributed around two values as shown in [19]. These two power law index 
values may be given a physical interpretation: the higher value (~ 2.1) stems from a 
shock- accelerated distribution of electrons and the lower value (~ 1.5) from the cooling 
electrons originating from the same distribution. The broken power-law component is 
described by: 

N(Et) -/^(^W"(t^)''' ^<^^W 
^ - \A{t),,,{l^r-P^ (^)-P^ E > E,{t) 

A high-energy cut-off was assumed due to the competition between the acceleration of 
the electrons and the radiative cooling that leads to a maximal energy that the electrons 
can be accelerated to. The values chosen for the simulations were based on results by 
de Jager et al. [20]. Their studies show that the spectrum of the Crab nebula has a high- 
energy cut-off e-folding energy of ~30 MeV, as predicted by their model. This is a 
robust value given by the model and independent of the magnetic field strength, B. We 
assumed, in our hybrid model implementation, that similar physical conditions exist at 
the GRB radiation site. Since the spectra from the prompt emission is boosted with a 
Lorentz-factor F ~ 100, we selected e-folding values of 5.0 GeV for GRB911016 and 
3.0 GeV for GRB941026 and GRB960530. These values fall within the GLAST energy 
range and the extended hybrid model is therefore in the simulations described by: 

^(P ,.\Nbb{E,t)+Nbpo{E,t), E<Ec 
^ " \Nbb{E,t)+Nbpo{Ej)e^'^^-^)l^^, E>E, 

where Ec is the energy from where the high-energy cutoff starts and E^ the e-folding 
energy. Extra- galactic background light (EBL) is expected to contribute to a high- 
energy cut-off for GRBs, but we did not consider this 77 absorption process in our 
simulations. This absorption process may however be important in spectra from GRBs 
at large cosmical distances. 
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TABLE 2. Temperature parameters used as input to SBM 







GRB911016 


GRB941026 


GRB960530 




[keV] 


67.23 


67.6 


37.38 


tn 


[s] 


0.992 


0.960 


2.392 


to 


[s] 


1.24 ±0.06 


2.1 ±0.2 


3.0±0.3 


a 




-0.04 ±0.03 


0.0 


-0.22 ±0.05 


b 




-0.48 ±0.02 


-0.35 ±0.02 


-0.65 ±0.04 


5 




0.08 


0.2±0.1 


0.02 



2.2. Modeling the temporal evolution 

Kocevski et al. [21] show that the FRED light-curves associated with the temporal 
evolution of the flux of single pulse bursts may be described by: 



F{t)=Fj ^ 



m 



d r f t 



+ 



d-\-r d + r\t, 



r+l- 



-{r+d)/{r+l) 



(6) 



where is maximum flux, the time at maximum flux, r the dimensionless power-law 
index describing the rise phase and d the corresponding index for the decay phase. The 
power-law describing the rise phase is in equation (6) proportional to f while the decay 
phase is described by: 

as suggested by Ryde and Svensson [22, 23] for single pulse bursts. Fq is here the 
flux at the start of the decay phase and T a normalization parameter. By utilizing the 
Levenberg-Marquardt algorithm we extracted the r and d parameters using equation (6) 
for the BATSE flux of each model component (blackbody and single power-law) over 
the energy range 20 keV-1 MeV for the three bursts. 

Ryde shows in [12] that kT in equation (3) is a function of time that often follows a 
smoothly broken power-law with an initial power-law index a and a later index b: 

\tnj \COSh[\og{tn/tQ)/5]J 

where to is the time of the break, 5 is the width of the transition, ^ = {b — a)/2, and 
^ — {b+a)/2. Equation (8) was used to extract the temporal parameters associated with 
the blackbody component in the hybrid model. The results from the fitting procedure, 
also here using a Levenberg-Marquardt method, are presented in table 2. 

The power-law indices, for all three bursts, exhibited some form of temporal evo- 
lution, but it was - due to the large variance for each individual data point - hard to 
determine whether the index was smoothly changing over time or in discrete steps. One 
interpretation of the varying power-law index is that the break in synchrotron radia- 
tion spectrum - that stem from a shock-accelerated electron distribution (higher energy 
power-law index) and from cooling electrons (lower energy power-law index) - moves 
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or jumps across the BATSE energy window as the burst progresses. For GRB911016 
and GRB941026 we simulated the temporal evolution of the power-law break with the 
following equation: 

i^n.O ^n.oo ) ^ ^ " ~\~ Efi aa (9) 

where Enfi is the initial energy at time t = Q, En.oo the asymptotic final energy when t 
tends to infinity and f„ the rate scaling factor. This equation is neither empirically nor 
analytically deduced, and is only motivated by it being a candidate for describing the 
observed evolution. The power-law component in GRB960530 could be well fitted with 
a constant value and we did therefore not implement any temporal evolution of its power- 
law break. We implemented, however, for this simulation of GRB960530 a fixed break 
at 6.0 MeV, with a high-energy power-law index of 2.1, as shown in figure 2. 

3. SIMULATING THE THERMAL GRBS 

A generic C-i-i- framework. Simple Burst Modeler (SBM), was developed in order to 
simulate the spectral evolution of the bursts covering the GLAST energy range. SBM 
produces histogram files, one for each detector-type on-board GLAST, that describe the 
amount of photons produced by the simulated burst binned in energy and time. The 
right panel in figure 1 show the time-resolved spectrum of the simulated GRB9 11016 
while figure 2 shows burst data simulated for GRB960530. The histogram file produced 




FIGURE 2. Evolution of the v/v spectrum over time in the keV-GeV range for a simulated GRB. 
The simulation is based on sub-MeV data collected by BATSE onboard CGRO for GRB960530. The 
blackbody peak can be seen around 100 keV, the power-law break around 10 MeV and the high-energy 
cutoff in the GeV domain. 

for the LAT instrument is then fed into gtobssim, the LAT fast observation simulator, 
and the Nal and BGO histogram files into GBM Tools, the GBM simulator. A modified 
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version of gtobssim was used so that generic histogram files of the type produced by 
SBM could be read in and used in the simulations. When running gtobssim we used the 
LAT instrument response file produced for Data Challenge 2 (DC2) in our simulations. 
GBM Tools was also modified with the added capabilities to output the detector energy 
grids into files required by SBM and to read in the histogram files produced by SBM. 
Figure 3 summarizes the steps that were performed for each simulated burst. These 



BATSE Trigger Data 



GLAST Definition fiie 



XSPEC 



Model Parameter Values + guesses 



GBM Simulator 

I 



Nal/BGO Energy Grid 




SBM 




LAT lightcurve histogram file 



Separate Nai and BGO lightcurve histogram files 



gtobssim 

gtselect 
gtrspgen 



i 



LAT FITS and response fll 




gtbin 



GBM Simulator 

I 

Nal/BGO TTE, background and response files 



® 



XSPEC 



® 



FIGURE 3. Combined SBM and GLAST simulation procedure. Note that gtobssim has to be executed 
before the GBM Simulator since the definition file required by the GBM simulation package is created by 
the GLAST ScienceTools in connection to the LAT simulation run. 

seemingly awkward steps were needed since gtobssim produces a definition file used as 
input to GBM tools and GBM tools create energy grid files required by SBM. 

1. XSPEC was first used to analyze the BATSE data in order to extract the necessary 
parameters used as input to SBM. 

2. SBM was executed to produce a LAT histogram file. 

3. gtobssim used the LAT histogram as input to produce the LAT FITS and response 
files, as well as the burst definition file, gtobssim also needs a response function and 
a template that describes the position of the burst relative to the GLAST spacecraft. 
During our simulations the response function for DC2 was used. We also utilized 
the same template in all simulations hence positioning all three bursts at = 62.9° 
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and = 95.3°, where and represents the inclination angle to the normal of 
the LAT detector and azimuthal angle around the LAT normal respectively, gtselect 
was used to select the events of interest and gtrspgen to produce the response files 
readable by XSPEC. 

4. gtbin was executed in order to create the FITS files in the PHA format readable by 
XSPEC. At this point all LAT data was produced. 

5. The GBM Simulator was now executed with the GLAST definition file as input, at 
this time only producing Nal and BGO energy grid files. 

6. SBM was executed again to produce histogram files that was used as input to the 
GBM simulator. The same parameter data was used as previously for the LAT run, 
but this time with the additional BGO and NaT energy grid files as input. 

7. The GBM Simulator was fed with the SBM histogram files producing Time- 
Triggered Event (TTE) files as well as background and response files. 

8. gtbin was now utilized again but this time on the files produced by the GBM 
simulator converting them to the PHA format readable by XSPEC. 

9. The produced LAT data and GBM data were now jointly analyzed with XSPEC. 

4. RESULTS 

XSPEC was utilized on the data produced by gtobssim and GBM tools by imposing 
the extended hybrid model, as described by equation (5), on the time-integrated spectra 
from two NaT detectors, one BGO detector and the LAT instrument. Both and C 
statistic, a modified version of Cash statistic [see 24, 17], was used in our analysis 
to estimate parameters that covered the first five seconds from trigger, while only 
statistic was used for the model test. C statistic give better results than X'^ at count 
rates below 10 counts per bin, but can not, as opposed to x^ tests, be used to get a 
"goodness-of-fit". It is therefore primarily used for estimation of parameter values. C 
statistic also assumes that the error on the counts is pure Poissonian, and should hence 
be the preferred statistic when estimating parameter values based on data with low 
photon counts. At higher count rates, both x^ and C statistic are expected to give the 
same parameter estimates, assuming we are dealing with Gaussian distributions. The 
model analysis was performed on both ungrouped and grouped data. In the latter case, 
the data was grouped into bins with at least 10 photons per bin for all detectors. Some of 
the results, from our analysis using x^ and C statistic, are presented in table 3. Figures 
4 and 5 show a simulated photon spectrum, based on BATSE data for GRB911016, 
together with the fit of the extended hybrid model. The left panel in figure 6 shows the 
vFy spectrum for the same burst. The resulting fit parameters were, for each simulated 
burst, compared with the SBM input parameters. The simulated GRB911016 output 
parameters was consistent with the SBM input parameter values. Table 3 presents some 
of the parameter values and also the reduced X'^- (v represents the degrees of freedom) 
and P- values from the hybrid model fits. All the fits showed reduced x^ values around 
1.0 and with high P- values, close to 1.0, for the ungrouped data. The high P- values 
indicate that the errors for the ungrouped data points, used in our XSPEC fits, are over- 
estimated. Grouping the data into at least 10 photons per bin, decreased the P-values 
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TABLE 3. Parameter values from hybrid model fit on simulated GLAST data 



Burst 


Statistic* 


Blackbody 


Powerlaw 


Powerlaw 


Model Statistics 






Temp [keV] 


Index 1 


Index 2 


V 


Av 


r 


GRB911016 




53.3± 1.1 


1 30 ± 04 


1 71 ±0 03 


336 


0.90 


0.92 






53.8± 1.1 


1 25 ± 04 


1 76 ± 03 


206 


0.98 


0.55 




c 


53.9±1.1 


1.26 ±0.04 


1.74 ±0.03 








GRB941026 




59.2±1.3 


1.6 (frozen) 


2.1 (frozen) 


339 


0.82 


0.99 




Xg 


61.5 ±1.3 


1.71 ±0.03 


2.9±0.8 


228 


0.82 


0.98 




c 


61.0±1.3 


1.69 ±0.02 


2.4 ±0.2 








GRB960530 


Xu^ 


41.8±0.8 


1.48 ±0.02 


2.4 ±0.2 


336 


0.91 


0.87 




Xu^ 


41.8 ±0.8 


1.42 (frozen) 


2.1 (frozen) 


339 


0.96 


0.70 




Xg 


41.7 ±0.9 


1.44 ±0.02 


2.3±0.1 


208 


1.03 


0.36 




c 


41.5±0.9 


1.42±0.01 


2.2 ±0.1 









* Xu indicates that model fit and parameter estimation was performed on ungrouped data, Xg that the 
same procedure was performed on grouped data with at least 10 photons per bin while C indicates that C 
statistic was used to estimate the parameter values. 



and also gave ;f^-statistic parameter estimates closer to the values found when using 
C-statistic. The right panel in figure 6 shows the confidence regions (with la, 2o 
and 3(7 levels) for the power-law photon-index versus the blackbody temperature kT. It 
shows a consistency between the BATSE data and the simulated GLAST data for the 
modeled temperature. The temperature was also consistent with the input data for the 
other two bursts, but the mean value for the photon indices became slightly softer in the 
GLAST data fits. This is probably due to the evolution of the blackbody temperature 
and, in the case of GRB941026, the evolution of the break-energy of the power-law. 
Freezing the indices to the expected values, when analyzing the ungrouped data, we still 
received good statistics ■^■> consistency with the input parameters. This can be 

seen for GRB960530 in table 3. It was however not possible to determine the confidence 
levels for the high-energy cut-off parameters in the fit - neither with binned x^ nor C 
statistic - due to the low photon count at the higher energies. 

5. DISCUSSION 

The hybrid model, that may be interpreted as a combination of photospheric and 
optically-thin synchrotron radiation, fits the spectra of single-pulsed GRBs well [see 
13, 19]. The problem with the synchrotron "line of death" (a > —2/3) is avoided with 
this model, since the Rayleigh- Jeans portion of the blackbody component explains the 
hard spectral slopes seen in many BATSE bursts. We presented here simulations of ther- 
mal bursts with a large super-MeV emission in the GLAST energy range using gtobssim 
and GBM tools. The simulations were followed by XSPEC analysis of the first five sec- 
onds of the resulting data for each burst, grouped to at least 10 photons per bin. This 
analysis showed that an applied thermal model, consisting of a blackbody function and 
a broken power-law with a high-energy cutoff, gave the desired "goodness-of-fit" repre- 
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GBM Simulation 




le-07 I ' ' ' ' ' ' 

10 100 1000 10000 

Energy [keV] 



FIGURE 4. Unfolded photon-spectrum from one Nal and one BGO detector of the GBM instrument 
covering the first five seconds of the simulated burst. The fitted two components of the hybrid model are 
also shown in this figure. The fit was performed on ungrouped data using statistic. The simulation was 
based on BATSE data for GRB911016. Note that the data points have been grouped in the plot. 



sented by the P-value close to 0.5 for all bursts. Parameter estimates of the power-law 
indices generally gave softer values than used for the simulation, but this can be ex- 
plained by the temporal evolution of the blackbody peak and power-law break over the 
five seconds that were integrated into one spectrum. The simulations of GLAST data 
for thermal bursts hence show that a photospheric component should be clearly detected 
by the GLAST instruments, if it dominates in the energy window of the GBM instru- 
ment and is super-positioned over an optically-thin synchrotron "background" spectrum 
extending from the lowest GBM energies into the LAT domain. The evolution of the 
Planck function will be possible to determine using events collected by the two detector 
types, Nal (8 keV - 1 MeV) and BGO (150 keV - 30 MeV), that constitute the GBM 
instrument. The non-thermal component will be detectable by both GBM and LAT. 

From the simulated GLAST data of the three single-pulse bursts it was not possible 
to determine the la confidence region for the high-energy cut-off parameters due to the 
low photon count, not even when using C statistic. We therefore expect that the cut-off 
and e-folding energy for some bursts may be hard to determine, especially for time- 
resolved spectra, due to the low photon count, if the LAT instrument response function 
is similar to the one used for Data Challenge 2 (DC2), even if the cut-off lies within the 
LAT domain. Possibly better results could be reached using the raw event data or using 
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le-12 



10 100 1000 10000 100000 le+06 le+07 le+08 

Energy [keV] 

FIGURE 5. The same simulated burst as in figure 4, including LAT data. The fit was performed on 
ungrouped data using statistic. Note that the data points has been grouped in the plot. 



a response function with more lean cuts than the ones used in DC2 for time-integrated 
spectra. 
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